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We calculate the particle production rate in an expanding universe with a three-torus topology. We discuss 
also the complete evolution of the size of such a universe. The energy density of particles created through the 
nonzero modes is computed for selected masses. The unique contribution of the zero mode and its properties 
are also analyzed. 



I. INTRODUCTION 



II. THREE-TORUS METRIC 



Although current astrophysical observations provide pre- 
cise information on the geometry of the universe 1 1 1, its topol- 
ogy remains a mystery. We don't even know whether the uni- 
verse is compact or infinite. Nevertheless, lower bounds can 
be put on its size for each compact topology (see, ll) and ref- 
erences therein). 

Amongst the possible topologies for the universe those with 
some or all spatial dimensions compactified are especially in- 
teresting, since then Casimir energies provide an additional 
contribution to the energy density. This may lead to an in- 
teresting vacuum structure for the standard model coupled to 
gravity which is insensitive to quantum gravity effects |3jt5j|. 
The simplest flat topology with all dimensions compactified is 
a three-torus, and that is the topology we will concentrate on 
throughout this paper. 

The scenario of our uruverse having a three-torus topol- 
ogy was investigated by many authors (see, |6l and references 
therein). Probably one of the most appealing features of such 
a model is that the creation of a three-torus universe is much 
more likely to occur than that of an infinite flat or closed uni- 
verse |6, 7|. In addition, it has also been shown that a three- 
torus topology can provide convenient initial conditions for 
inflation [7J. 

Here we consider gravitational particle creation in an ex- 
panding toroidal universe. The particle production formalism 
was developed in liSnlOl and investigated in great detail in 
later works (see, lfTTI - [T3l and references therein). Since then, 
it has been thoroughly studied in the case of a FRW cosmol- 
ogy, including its implications for dark matter creation around 
the inflationary epoch lfT4ll . However, particle production in a 
toroidal universe hasn't been extensively studied (see, (TSl for 
some work on the subject). 

In this paper we provide a detailed numerical calculation 
of the particle production in a three-torus universe. We start 
with introducing the relevant formalism. We then discuss the 
evolution of the size and energy density of the universe from 
the Planck time to the present time. Next, we derive analyt- 
ical formulae for the particle number and energy density at 
early and late times. We then find full numerical solutions and 
confirm that they agree with the analytical approximations in 
the appropriate regions. The particle production through the 
nonzero modes is somewhat similar as in the closed universe 
case discussed in 1 12|. However, the three-torus particle cre- 
ation includes an additional contribution from the zero mode, 
which strongly depends on the choice of initial conditions. 



We start with the following spacetime interval, 

ds^ = -dt^ + Ujdy'dy' , (1) 

where tij is the metric on the three-torus with i,j = 1,2,3, 
and the compact coordinates are ?/* G [0, 2tt). The 3x3 matrix 
with components tij is positive definite and has a determinant 
equal to the volume modulus a^. A suitable parametrization 
is given by. 




Pi 

Pin + P2T2 
pI + pI + pI 



(2) 



where (ti, T2, pi, p2, Pa) are the shape moduli. We assume 
that all the parameters in (|2]l are independent of the spatial co- 
ordinates. Furthermore, we seek stable solutions to Einstein's 
equations for which the shape parameters are also constant in 
time. It was shown in |5] that this is only possible for. 



{n,T2,Pl,P2,P3) 



2' ~Y' 2' ~6~' ~Y 



(3) 




which arises from the symmetries of the Casimir energies. For 
a three-torus universe characterized by the parameters ([3]) only 
the evolution of the volume modulus is nontrivial. The corre- 
sponding metric takes the form. 



(4) 



and this is the metric we will adopt in our further analysis. 



III. PARTICLE PRODUCTION 

We first derive general formulae for the gravitational parti- 
cle production rate in a three-torus universe. Let us consider 
a complex scalar field = 'i'{x) of mass m with the La- 
grangian density given by. 



(5) 
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where R is the Ricci scalar, 
6 



R = ^ (a d + d'^ 



(6) 



The dot denotes the derivative with respect to time t and the 
factor ^ = 1/6 was chosen to have conformal invariance in the 
limit m 0. We note, however, that the results in this paper 
are not very sensitive to this choice and would be similar, for 
example, in the case of a minimally coupled scalar field (for 
which $ = 0). 

The stress-energy tensor is given by. 



1 



c 



The equation of motion for the field 4' is. 



a 



a 



f = 



(7) 



(8) 



where is the Laplacian on the three-torus. We write the 
solutions of equation ([8]) as. 



with A = (Ai, A2, A3). In our case, 

i,. 



AiAj0A, 



where. 



(9) 



(10) 



(11) 



and there is an implicit sum over the repeated indices i,j ~ 
1, 2, 3, with Ai = 0, ±1, ±2, .... The formula for 0^ is given 

by. 



(12) 



Note that for A = (0, 0, 0) we have (po — const, which corre- 
sponds to the zero mode. Equation ((8]l takes the form. 



a \ a 



a 



ux = , 



where. 



(ma)^ + iijXiXj 



(13) 



(14) 



We can now quantize the field introducing standard canonical 
equal-time commutation relations for the field and its general- 
ized momentum. Those relations are satisfied if we write the 
field as. 



* = 



1 



(2^a)V2 



^ (pX ul OA + (j)*x ux b{ 



(15) 



Ai ,A2 ,A3 — — 00 



where d\ and dx are the creation and annihilation operators of 

a particle in the state A = (Ai , A2 , A3), and b\ and 6a are those 
for antiparticles, all obeying the usual commutation relations. 
Adopting such a convention, the Hamiltonian can be written 
as I12i, 



Ai,A2,A3 — — 00 



WA 
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B*xdxbx 



(16) 



with. 
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Bx 



2ojx 
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(17) 



and A defined through (px — 



In general, the Hamiltonian ( I61 is non-diagonal. The re- 



quirement of it being diagonal at some initial time <o imposes 
the following conditions. 



uxih) 



1 



\/wA(to) 



uxih) 



\/^A(to) 

a(io) 



(18) 



Now, we can diagonalize the Hamiltonian through the follow- 
ing BogoUubov transformation, 

dx = al{t)d\{t)+Pxit)b'lit) , 

bx = a*x{t)b'x{t) + ^xit)alit) , (19) 



where |q;aP ^ l/^Ap = 1- It can be shown lfT2l . from the re- 
quirement that there be no non-diagonal terms d'^b'^ or d'J)'s^, 
that the equations for ax{t) and f3x{t) are. 



« '^A / „. f'ujxit', 

Pa = CkAexpl— 2i / — r-r^dt 



ax 



2WA V ito «(^) 

/^Aexp 2i / 



(20) 
(21) 



with the initial conditions /3a (^o) and axito) — 1. The 
function ux is expressed in terms of a a and /3a as [121, 



MA = 



QfA exp I I I — 7-7:- dt 



to 



a{t') 



■ /3a exp 



a{t') 



dt' 



(22) 



Relations ( |20] i and ( |2T| ) can be combined into one second order 
differential equation for /3x{t), 



/3a + 



WA 



t^A 



2iujx 



^ujx ojx a 
with the initial conditions, 

/3A(io)=0, /3a (io) 



/3a 



/3a - , (23) 
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(24) 
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epoch 


pre-inflationary era 


inflation 


radiation era 


matter era 


dark energy era 




^ 5 • 10-*" s 
~ 8 - lO-^'^ GeV-^ 


^5 - 10-^*^8 
~ 7 - 10-12 GeV-i 


^ 10"^=* s 
~ 10-" GeV-* 


^7-10* years 
~ 3 - lO^*' GeV-i 


^ 10*" years 
~ 5 • 10** GeV-* 


t'hnal 


^ 5 ■ 10"^" s 
~ 7 - 10-12 GeV-i 


10"^^ s 
~ 10"'' GeV-i 


~ 7 • 10* years 
~ 3 • 10^^ GeV"* 


^ 10*" years 
~ 5 X 10** GeV"* 


^ 1.37 • 10*" years 
~ 7 • 10** GeV"* 


a{t) 


, ^ „ 1/1 1 






a„o ~ 6 - 10'*^;^ 




ttio exp [^ptot/(3M2) 
a,o ~ 10"!^ GeV"* 


fldo oxp [^Ptot/ (3M2) t] 
Qdo ^ 2 • 10*2 GeV-* 






V GcV 


a(tiiiit) 


2 • 10"^^ m 
~ IQ-i^ GeV-i 


^ 2 ■ IQ-^'i m 
~ 10-1^ GeV-i 


^ 4 m 
~ 2 • 10*^ GeV-i 


« 2 ■ 102^ m 
~ 10^" GeV-* 


8 • 102'^ m 
~ 4 • 10*2 GeV"* 


a(^flnal) 


2 • 10-^1 m 
~ 10"^^ GeV-i 


« 4 m 
~ 2 ■ 10*'' GeV"* 


« 2 ■ 102^ m 
~ 10^^ GeV"* 


« 8 ■ 102^ m 
~ 4 - 10*2 GeV"* 


^ 102^ m 
~ 5 • 10*2 GeV-i 


Ptot(i) 




const 


3M2/(4t2) 


4M2/(3t2) 


const 


Plot l,tinit ) 


^ 7 • lO'^* GeV* 


^ 9 ■ 10^** GeV* 


^ 4 • 10"* GeV* 


^ 8 • 10-^^ GeV* 


« 3 • 10-*^ GeV* 


Ptot(tfinal) 


^ 9 • 10=** GeV* 


^ 9 ■ 10^** GeV* 
after reheating: 
^ 4 • lO"^* GeV* 


^ 5 • 10"^^ GeV* 


^ 3 • 10"*^ GeV* 


^ 3 • 10"*^ GeV* 



TABLE I: Estimated timeline for the evolution of a three-toras universe from the Planck time until the present time tuniv ~ 13.7 billion years, including the 
size of the universe and total energy density during different epochs, assuming a pre-inflationary expansion a(i) = apo\/i and the current size of the universe 
tt(tuniv) ^IORh- 



The normal ordered Hamiltonian now takes the diagonal form, 

oo 

H= ^x{all^'x + yiK). (25) 

and the physical vacuum depends on time through, 

OaWIOW) -6^(t)|0(i)> =0. (26) 

It is now straightforward to arrive at the formula for the num- 
ber density of particle pairs created after time t. 



n{t) 



1 

(27ra)3 



(27) 



The energy density and pressure of the created particles are 
calculated from the vacuum expectation values of the appro- 
priate components of the stress-energy tensor (|7|. After nor- 
mal ordering the result is (compare with [1121 '). 



Pit) 
Pii{t) 



(28) 



Al ,A2 ,A3 — — OO 



127r3a4 
Pii{t) 



E 



Ai,A2,A3 — — oo 



for i ^ j , 



(29) 
(30) 



where i.j = 1,2,3. Note that the pressure is not isotropic in 
our case. It would be isotropic only for the simplest choice of 
the three-torus metric diag(a^, a^,a^). 



IV. EVOLUTION OF THE THREE-TORUS UNIVERSE 

In order to calculate the particle production rate, we first 
have to know how a{t) evolved in time. Our case is different 
from the usual FRW cosmology, in which a{t) is the scale 
factor and can be rescaled by an arbitrary number For a three- 
torus universe a{t) has a physical meaning - it describes the 
size of the universe at a given time t. Recent analyzes (see, 
|fT6l and references therein) set a lower bound on its present 
value of a(iuniv) ^ QRh, where Rh is the Hubble radius 
today. Let's therefore assume that the current size of the three- 
torus universe is. 



a(iuniv) = 10 Rh 



10^^ m 



(31) 



Table I shows the estimated timeline, the size of the uni- 
verse and the total energy density during different epochs in 
the evolution of the universe. For each epoch the total energy 
density is uniquely determined by the expansion rate through 
the Friedmann equation. 



Ptot 

3M2 



(32) 



where the reduced Planck mass Mp ~ 2.4 • 10*® GeV. Us- 
ing the formulae for a{t) for different epochs, we can evolve 
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a(iuniv) back in time and find its value at any particular in- 
stance in the past. 

Observational data suggest that we currently live in a 
dark energy dominated universe, where the energy density is 
roughly constant and the expansion is exponential. It was pre- 
ceded (at t < 10 billion years) by a matter dominated era with 
the size of the universe increasing like t^^^. Before the matter 
era (at t < 70000 years), radiation was driving the expansion 
like Vi. It is believed that before the radiation epoch there 
was a brief period of inflation (lO"^'"^ s < i < lO^'^"^ s), an 
exponential expansion resulting from a constant energy den- 
sity. The inflationary energy density in table I was estimated 
adopting 72 e-folds of inflation and assuming that it ended at 

t W 10^33 

The expansion rate in the pre-inflationary era is unknown. 
However, in the case of a three-torus universe it might have 
a natural explanation. As mentioned in the introduction, a 
compact topology results in the appearance of Casimir ener- 
gies of existing fields. The full formulae for the Casimir en- 
ergies in a three-torus universe for an arbitrary set of shape 
moduli was derived in |5|. For a real scalar field with mass 
m <C 10^^ GeV the Casimir energy density before inflation 
for our choice of metric is. 



(27r)4 a4 I 22/3 12 ^ 32/2 tt 360 

00 ^ s 3/2 

33/4 X! \^) COS (tt n2 na) K3/2{y^ 7^2 713^ 



PCas(a) 

211/6 



"2. "3 = 1 



+ ^ COS Ittoi (712+^3) Y ^y^(n|-n2ri3 + n§) 

"2,"3=— 00 ni — 1 



^ Ki (^^7rniy^(^^2T^^^^^^r,r^^ 
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FIG. 1: Size of the three-torus universe as a function of time, assuming 
a(tuniv) = 10 Rh- 
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FIG. 2: Total energy density in the three-torus universe as a function of time. 

Started, 



which can be written as. 



yOCas(a) 



K 



(34) 



where k~5x10^'^. In a general case, this formula in- 
cludes also a factor corresponding to the number of degrees 
of freedom for a given field and an additional minus sign for 
fermions. Interestingly, the Casimir energies at a <C 1/m 
satisfy the same equation of state as radiation, i.e.. 



PCas(a(i)) = ^pcas(a(t)) • 



(35) 



Therefore, if there were more fermionic degrees of freedom 
{nf) than bosonic (n?,) and Casimir energies dominated the 
total energy density before inflation, the universe would ex- 
pand according to a(i) — ao\/i, with the total energy density 
given by. 



tot 
PCas 



(«) 



{nf - nb)K 



(36) 



Surprisingly, if we choose in our case {nj — rif,) ~ 20, we ob- 
tain the total Casimir energy density at the time when inflation 



PCas(a(^inf)) «Ptot(iinf) 



(37) 



where a(tinf ) is the universe size at the beginning of inflation 
and ptot (^inf ) is the total energy density of the universe during 
inflation, both of which where estimated before independently 
of the Casimir energies through the evolution back in time. Of 
course, the required value of (n/ — rib) which satisfies condi- 
tion (37 1 strongly depends on the inflationary parameters. 

Another interesting observation concerning the three-torus 
topology is that assuming the current size of the universe 
a(^univ) ~ 10i?^f and a pre-inflationary expansion a{t) = 
dpoVi, the evolution back in time yields the size at Planck 
time a{tp) « Ip, where Ip is the Planck length. Note that an 
evolution from the Planck size at Planck time to ten Hubble 
radii at present time wouldn't be possible if we assumed a 
universe with a closed topology, since then the curvature term 
contribution would entirely dominate the total energy density 
during the early stages of the universe evolution. 

The plot of the size of the universe a{t) corresponding to 
the values from table I is given in figure 1 . A similar plot for 
the total energy density is shown in figure 2. We note that a 
similar calculation to the one presented in this paper can be 
done assuming a different pre-inflationary expansion. 



5 



Knowing the shape of a{t), we can now numerically solve 
equation ([23]l for the Bogoliubov coefficients and then use for- 



mulae ([27[i and (28 1 to calculate the number density and en- 



ergy density for the gravitationally created pairs of scalar par- 
ticles of mass m. Our results can be easily generalized for 
other types of particles by adopting the appropriate form of 
the stress-energy tensor We will perform the calculation for 
three different masses: mi = 10^ GeV, 1112 = lO'^ GeV, and 
TO3 — 10~^ GeV, without the zero mode first, and then com- 
puting its contribution as well. In our calculation we will be 
using natural units, i.e.. 



2 • 10 



-16 , 



1 GeV" 



7-10 



-25 



(38) 



We will also assume that there is no back-reaction of the cre- 
ated particles on the background evolution. 

V. ANALYTICAL RESULTS FOR NONZERO MODES 

In order to fully understand the numerical solutions it is 
very helpful to derive their analytical behavior in two regions: 
t <C 1/m and t ^ 1/ni. In contrast to [11, 12), we choose 
not to specify the initial conditions at the singularity, but at the 
Planck time, since the physics at earlier times is unknown. 



(39) 



However, we note that for the nonzero modes the solutions 
are insensitive to the value of to at which we impose those 
conditions, as long as we look at the region t ^ Iq. 

Let us first consider the case mt ^ 1 and focus only on the 
pre-inflationary epoch, for which we adopted a{t) — apoVt. 
With the additional assumption ma <C 1 the largest contribu- 
tion to the sum ( [27] i comes from terms with small A, . Thus, we 
can assume Xit <^ a, in which case the solution to equation 



(23 I is given by. 



(iijXiXj) 1^ 



(40) 



After plugging (HOll into ([27|, the particle number density is. 



n{t) 



1 



(27ra)3 



E' 



•^1 -'^2 5^3— ~ 00 



1287r3 



Ai,A 



00 

E 

.A3 = 



-00 i^iJ^i^lY 



250 



(41) 



where the prime excludes the zero mode. 

The calculation of the energy density and pressure in the re- 
gion mt ^ 1 cannot be performed with the same assumption 
Ait <C a, as now terms with large A; contribute significantly. 
However, for ma <C 1 we have \f3\\ ^ 1 and ~ 1. With 
equation ( 20 1 we can approximate I3\ (t) by. 



/3a 



2 tijXiXj jQ 



a{t')a{t') 



X cxp 



dt' 

a{t') 



(42) 



Using this result to rewrite ( 28 1, we arrive at the energy den- 
sity for the created particles during the pre-inflationary era. 



p{t) 



4 4 
647r3a4 



E ' 



1 



Ai ,A2,A3 — — 00 



-00 i^^i^i^jy^^ Jo 



dti / dt2 



X exp 



-2i{i,jX,XjY/'— (Vti - Vh) 



(43) 



The biggest contribution in equation ( |43] l comes from sum- 
ming the terms with large A^. It turns out that we can ap- 
proximate our three-torus by a three-sphere of radius \^a, 
for which it is possible to perform the sum over A's. Equation 
( |43| ) takes the form. 



Pit) 



4 4 



00 ^ 



exp 



A=l 



t pt 

dti / dt2 
Jo 



Giy^n^a^Jo Jo 



dti / dt2 log 



1 



1— exp 



apo \ V 



(44) 



The double-integral above is real and positive. The resulting 
energy density for mt ^ 1 is essentially constant in time. 

Using a similar approximation one can derive the formula 
for the pressure. 



Pijit) 



r Pit) 



t, 



4 2 
_ ma 



'pO 



y 3 

X Lis ( 6xp 



967r3a2 
" 2i 



t pt 
dti / dt2 



2Re dti log 
'0 



1 — exp 



2i 



iVh-Vt2) 
h-Vi 



(45) 



A numerical check shows that for mt ^ 1 the pressure for 
the created particles is also almost constant and satisfies the 
"quasi-vacuum-like" equation of state. 



Pijit) ~ pit) 



(46) 



We will not go into the details of calculating the parti- 
cle number or energy density in a three-torus universe for 
mt ^ 1, since this case is analogous to that of an infinite 
flat or closed universe, which was explored in detail in |12|. 
Briefly summarizing, for mt ^ 1 all Bogoliubov coefficients 
Px are roughly constant and particle creation doesn't occur. 
For a universe expanding according to a{t) = a^f, where 
< g < 2/3, the number density and the energy density of 
the created particles behave in the following way. 



n(t) 
Pit) 



t^l 



1 



1 

t39 ^ a^ 
while the pressure Pijit) <C pit). 



(47) 
(48) 
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FIG. 3: Number density of the created particles for mi = 10^ GeV (blue 
solid line). The green dashed line corresponds to the mt <C 1 and ma <g; 1 
approximation given by equation \4i\ and the orange dot-dashed line corre- 
sponds to the mt S> 1 approximation given by equation J47j. 
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FIG. 4: Energy density of the created particles for mi = 10^ GcV (blue 
solid line). The green dashed line corresponds to the constant density for 
mt <C 1 given by equation j44j and the orange dot-dashed line corresponds 
to the mt S> 1 approximation given by equation {41) . For comparison, the 
total energy density of the universe is plotted in solid red. 



VI. NUMERICAL RESULTS FOR NONZERO MODES 

Having an analytical understanding of how the functions 
n{t) and p{t) behave for mt <C 1 and mt ^ 1, in this section 
we discuss the full numerical solutions. 



We first solve numerically equation (23 i and find l3\{t) for 
different sets of Ai's, with bj\{t) defined through equation 
( 14 1. The initial conditions are chosen according to (39)1. We 



then plug the calculated /3a (t) to equations ( 27 1 and ( 28 



and 

in this way obtain the particle number and energy density. 

Figure 3 shows the plot of the particle number density 
(in blue) corresponding to a mass of the created particles of 
TOi = 10^ GeV. The green dashed line increasing like \/t 
is the mt <C 1 and ma <C 1 approximation given by equa- 
Since the particle mass is large, this approxima- 



tion (41 



tion breaks down relatively quickly and the number density 
starts slowly decreasing. During the inflationary epoch, it falls 
down rapidly by many orders of magnitude. Particles of mass 
10^ GeV are no longer produced after inflation. This is con- 
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FIG. 5: Same as figure 3, but for m2 = lO'^ GeV. 
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FIG. 6: Same as figure 4, but for m2 = 10^ GeV. 

firmed by fitting the line falling off like t^"^/^ (orange dot- 



dashed line) corresponding to the relation (47 1. 

Figure 4 shows the time evolution of the energy density for 
the created particles of mass TOi = 10^ GeV. The energy den- 
sity is nearly constant for mt ^ 1 and its value agrees with 



that from equation (44 1, represented by the green dashed line. 



The subsequent behavior of the energy density is similar as in 
the number density case - following a mild decrease there is 
a rapid fall during inflation, after which the density continues 
to decrease like ^^■^/■^ (orange dot-dashed line), since particle 
creation no longer takes place. Note that the density of the 
created particles is small compared to the total energy density 
in the universe (denoted in figure 2 by the red solid line) at 
any given time. 

Figures 5-8 show the number density and energy density 
for the created particles with masses TO2 — 10'^ GeV (figures 
5, 6) and m^ = 10~^ GeV (figures 7, 8). There are two qual- 
itative differences compared to the previous plots. Because 
the particles are lighter, the conditions mt <C 1 and ma <C 1 
hold also for later times and the plots match the approxima- 



tions (41 1 and (44i (green dashed line) all the way until the 
inflation epoch. Secondly, particle creation continues after in- 
flation and stops only at t 1/m, which is confirmed by 
fitting the line (orange dot-dashed) going like t~^^'^. 
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FIG. 7: Same as figure 3, but for r?i3 = 10^'^ GeV. 
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FIG. 8: Same as figure 4, but for ms = 10~^ GeV. 



We emphasize that although particle creation lasts arbitrar- 
ily long for light enough particles, their energy density is still 
negligible compared to the total energy density. For instance, 
if we wanted to have particles created today, we would need a 
particle of mass JTiiight ^ 10^*^ GcV. From formula (44 1 the 
energy density is suppressed approximately by m''/a^Tak- 
ing into account also the huge decrease in energy density of 
the created particles during inflation, there is no way it can 
ever compete with the total energy density of the universe. 



VII. ZERO MODE 

Let us now consider the contribution of the zero mode. In 
this case all the formulae simplify significantly. Relation ( 14 1 
becomes, 



UJQ{t) — m a{t) 



(49) 



Equation ( 23 i for the Bogoliubov coefficient reduces to. 



a a 



/3o + ( 7 - T + ) /3o - f^/3o = , (50) 



with the initial conditions, 

/3o(io) = 0, po{to) 
The equation for wo(i) is. 



ajto) 
2a(io) 



a / 2 o a 

uo + 3-uo + [ m H \ — - uo = 

a V a a'' 



where. 



uo{to) 



\/ma{to) 



a{to) 



(51) 



(52) 



(53) 



Once equations ( |50| l and p2} are solved, the particle num- 
ber density, energy density and pressure can be calculated us- 
ing the simple relations. 



l/3o| 



nit) 



(27ra)3 



Pi jit) 



(54) 
(55) 

ima|uoP + ^ j . (56) 



It turns out that equation ( |50] l with the initial conditions ( 5 1 
can be solved exactly for simple choices of ait). For instance 
in the pre-inflationary era with a(i) = ao^/i, the solution for 
a particle of mass m takes the form. 



U{ll2imto) 



X Li4%(2imt) - C/ (|, |, 2m7i) L^J^/^i2imto) 
4U{ll2im.to) L%%(2jmio) 



-5t/(|,|,2imto) L^J^/^i2tmto) 



(57) 



where I7(x, y, z) is the confluent hypergeometric function of 
the second kind, and L^ix) is the associated Laguerre poly- 
nomial. 

In the limit rnt <C 1 the dependence of the Bogoliubov 
coefficients in formula (57 1 on t and to is /3o ^ it/ta)^^*, 
which translates to. 



nit, to) pit, to) 



1 



(58) 



4a^ 



Thus, setting initial conditions at a small enough to, one can 
have an arbitrarily large zero mode particle production at any 
given instance of time. As mentioned earlier, this situation 
doesn't occur for the nonzero modes, for which $\ito) is in- 
dependent of to for ma ^ 1. However, since we don't know 
the physics before the Planck time, choosing to < tp is not 
justified. In the subsequent calculation we assume tg = tp, 
just as in the case of the nonzero modes. 
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FIG. 9: Number density of the particles created through the zero mode for 
m = 10'' GoV (blue solid line), assuming to = tp. The green dashed 
line corresponds to the mt <^ 1 approximation given by equation j58|. The 
orange dot-dashed line denotes the mt 2> 1 approximation given by equation 



> 

O 



o 




Log;o(t [GeV-']) 



FIG. 10: Energy density of the particles created through the zero mode for 
m = 10'^ GeV (blue solid line), assuming to = tp. The green dashed 
line corresponds to the mt <^ 1 approximation given by equation |58|. The 
orange dot-dashed line corresponds to the mt ^ 1 approximation given 
by equation |60l. The red solid line denotes the total energy density of the 



It turns out that also equation (52i possesses an analytical 
solution for a{t) — aoy/i, given by. 



Uo{t) 



Yi/iimt) 



4(mao\/t)i/2 
+ (1 + 2imio) Ji/4(mto) 



- 2mto J-3/4(™io) 
+ Ji/i{mt) 



2mtoY_3/4{m.to) - (1 + 2imto) Yi/4(mio) 



(59) 



where Ja {x) is the Bessel function of the first kind and (x) 
is the Neumann function. 

Figure 9 shows the number density and figure 10 presents 
the energy density for particles of mass ni — 10^ GeV cre- 



0.4 
0.2 
0.0 
-0.2 
-0.4 



p,,7p = - 2'/'/ 3 



p,,7p = 2'/3/ 3 



-15 -10 



Log,,,(t [GeV-i]) 



FIG. 11: Ratio of the pressure pii{t) and the energy density p{t) of the 
particles created through the zero mode for m = 10^ GeV assuming the 
initial condition at to = tp. 



ated through the zero mode, assuming a full evolution of a{t) 
given in table 1, with the initial conditions set at to — tp. In 
the region mt <C 1 both the number density and energy den- 



sity decrease according to (58 i, which is confirmed by fitting 



the line (green dashed) going like 1/t. Because of choosing 
a small mass of the created particles, the fit matches the plot 
until the beginning of inflation, during which both the num- 
ber density and energy density fall rapidly by many orders of 
magnitude. The orange dot-dashed Une corresponds to. 



n{t) ^ pit) 



1 



(60) 



i.e., no particle creation. It fits well the plot for t > 1/m, 
which implies that for the zero mode particle creation also 
stops at t ss 1/m. The energy density of the created parti- 
cles is again small compared to the total energy density in the 
universe. ' 

It is interesting to analyze the equation of state for the mat- 
ter created through the zero mode, as it has significantly dif- 
ferent properties than in the nonzero mode case. Let us for a 
moment assume that the entire evolution proceeds according 
to the pre-inflationary expansion a{t) = apQy/i. An analysis 
of formulae ( |57j l and (|59]l shows that for such a scenario. 



Prjit) = 



Uj \p{t) 



for mt <^ 1 
for mt ^ 1 



(61) 



Since ta — and t.^ = l/v^i for i ^ j, the created matter 
would have a "quasi-accelerating" equation of state for times 
t < 1/m. This picture, however, is significantly modified 



when we take into account inflation. 



We note, however, that if p{t) "C for a range of t < and the equa- 
tions in this paper are still valid in that region, by assuming a small enough 
to the energy density of the particles created through the zero mode might 
in general be large enough to compete with ordinary matter and radiation 
at some instance of time in the evolution of the universe. 
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Figure 1 1 shows the the ratio pa (t) / p{t) for the zero mode 
production of particles with mass m = 10'^ GeV including 
the full evolution of the three-torus universe from table I. As 
expected from equation ( [6T] ), for the pre-inflationary epoch it 
has the form pii{t) / p{t) = — ^^2/3 . However, it doesn't stay 
like this until t « 1/m - inflation itself changes the equation 
of state to pii{t) = p{t)^/2/Z, and it remains this way for 
the rest of the evolution. This behavior seems generic for all 
masses of the created particles. An interesting fact is that for 
much shorter inflation times one could end up with an equa- 
tion of state with Pii{t)/p{t) anywhere between -~\/2/3 and 
v/2/3 at the end of inflation. However, such short inflation 
periods are not realistic. 

Vni. CONCLUSIONS 

In the present work we calculated the rate for gravitational 
particle creation in a three-torus universe. We performed the 
calculation assuming the metric on the three-torus for which 
the shape moduli are stable. We adopted the current size of the 
universe of ten Hubble radii and estimated its full evolution in 
time. We argued that the Casimir energies might dominate the 
total energy density before inflation, resulting in the expansion 
of the universe like in the radiation era. We then showed that 
there are two types of contributions to the particle production 



rate - one from the nonzero modes on a torus, and the other 
coming from the zero mode. 

The nonzero mode contribution is interesting because at 
early times it is characterized by a "quasi-vacuum-like" equa- 
tion of state for the produced matter. The particle production 
itself continues until t Ri 1/m, after which it essentially stops. 
It turns out that the energy density of particles created through 
the nonzero modes is tiny compared to the total energy density 
in the universe at all stages of its evolution. 

The production through the zero mode is more unique. The 
corresponding energy density of the created particles is very 
sensitive to initial conditions. If we set those at the Planck 
time, the energy density of the particles is again small with 
respect to the total energy density of the universe. We showed 
that although the particles produced through the zero mode 
have a "quasi-accelerating" equation of state before inflation, 
they are described by a "quasi-radiative" equation afterwards. 
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